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Abstract 

As is well-known, there exists a four parameter family of local interac- 
tions in ID. We interpret these parameters as coupling constants of delta- 
type interactions which include different kinds of momentum dependent 
terms, and we determine all cases leading to many-body systems of dis- 
tinguishable particles which are exactly solvable by the coordinate Bethe 
Ansatz. We find two such families of systems, one with two independent 
coupling constants deforming the well-known delta interaction model to 
non-identical particles, and the other with a particular one-parameter com- 
bination of the delta- and (so-called) delta-prime interaction. We also find 
that the model of non-identical particles gives rise to a somewhat unusual 
solution of the Yang-Baxter relations. For the other model we write down 
explicit formulas for all eigenfunctions. 



1 Introduction 



By general physical arguments one expects that, to understand the low-energy 
properties of a quantum system, one can ignore details of the short distance struc- 
ture of the interactions and replace them by local interactions which are singular 
and non-trivial only at a point. The most prominent such local interaction is 
formally defined by a delta function potential and parameterized by one real pa- 
rameter, but it is known since quite some time that, in one dimension (ID), the 
most general local interaction is characterized by fourieal parameters While 
the earliest known examples of exactly solvable ID quantum many-body systems 
with two-body interactions are defined with the delta function potential |Sl IH E] 
(see also jE]), the full four parameter family of local interactions has received much 
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less attention in this context until recently [3 |HJ 13 CHI UH 1121 UHl HI] (see also 
Chapter 7 in ^^). 

In this paper we consider the ID quantum many-body systems with two-body 
interactions given by the general four parameter family of local interactions, and 
we determine all cases which are exactly solvable by the coordinate Bethe Ansatz 
even in the general case of distinguishable particles. (The same problem was 
also studied in [TU], but our approach and results are different, as discussed in 
more detail below.) We find and solve two families of such models which provide 
interesting generalizations of previously known cases. We use the parameterization 
of the general local interaction proposed in [12], and we suggest a natural physical 
interpretation of these parameters as coupling constants of delta-type interactions 
which also include different kinds of momentum dependent terms. This allows 
us to write down formal Hamiltonians to define and interpret these models in a 
simple manner. While this interpretation is different from others in the literature 
mE]; the mathematically precise formulation of our models is the usual one in 
terms of boundary conditions. In the following two paragraphs we describe in 
more detail the exactly solvable models which we find. 

The first family of models can be formally defined by the following Hamilto- 
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1r]b{xj - Xk) i {d^^ - d^^) (1) 

which depends on two real coupling constants c and t]. For this and all other 
many-body models considered in this paper we assume that the particles move 
on the full real line, Xj G M for j = 1, 2, . . . , A^, and we ignore bound states, 
which makes our solution complete only under certain restrictions on the coupling 
constants, such as c > 0.^ This model is a one-parameter extension of the famous 
ID delta gas solved, in the boson case, by Lieb and Liniger [3], and in the general 
case of distinguishable particles by Yang using the coordinate Bethe Ansatz [S]. 
We interpret the extension as a particular momentum dependent interaction.^ We 
find that it is possible to generalize Yang's solution to the full model, even though 
it describes identical particles only in the Yang case t] = (since the interaction 
terms with coupling r] are not invariant under particle exchanges). 
A formal definition of the second model is 

^ r 2 1 

H^"^ = - E + E i2c5(x, - X,) + - (a,, - a. J S{x, - Xk) {d., - a. J ] (2) 

i^i i<fc 

""^Bound states are of course interesting but beyond the scope of this work. 
^Since pj = —idx^ are the particle momentum operators. 
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with only one real coupling parameter c. The first interaction term is identical 
to that of the ID delta gas, but there is an additional momentum dependent 
interaction term with inverse coupling strength. As will be discussed below, this 
second interaction is identical with what is usually referred to as delta-prime 
interaction |21IIH|- This model describes identical particles, and it is thus possible 
to restrict it to bosons and fermions. It is interesting to note that, for fermions, 
the delta interaction is invisible (Pauli principle), and the model reduces to one 
first solved by Cheon and Shigehara [H], whereas for bosons the second interaction 
is invisible reducing it to the ID Bose gas solved in j3] (see also [12] )• We find 
that the explicit eigenfunctions of this model can be constructed in a remarkably 
simple manner by taking the well-known eigenfunctions of the boson delta gas 
restricted to the fundamental wedge Xi < X2 < ■ ■ ■ < and extending it to all 
other wedges xqi < xq2 < . . . < xqn, Q E S^, using whatever particle statistics 
one considers. This generalizes the duality found in ^ to arbitrary exchange 
statistics. 

We now recall some basic facts about general local interactions in ID (see e.g. 
|16j for a more comprehensive discussion). The simplest example of a model with 
such interactions can be formally defined by a Hamiltonian 

H = -dl + V (3) 

with X G M the particle coordinate and V an interaction localized at x = 0, 
i.e., the action on wave functions ip vanishes except at the origin. This implies 
that the Schrodinger equation determining the eigenstates ip of H is trivial nearly 
everywhere, ip" + Eip = for x ^ 0, but the interaction results in non-trivial 
boundary conditions at the singular point x = 0. Models of this type have been 
studied extensively using the theory of defect indices, from which it is known that 
the most general such interaction can be parameterized by four real parameters 
mi2]. The most prominent example is the delta function potential, 

V = c6{x), (4) 

parameterized by one real coupling constant c. Another well known special case 
is what often has been referred to as delta-prime interaction; see Section 1.4 of 
Recently it was shown that the boundary conditions defining this interaction 
arise from the following momentum dependent potential, 

V = AXdJ{x)d, (5) 

with A a real coupling constant y^.^ In this paper we extend this point of view 
to the full four parameter family of local interactions, and we show that they arise 

■^This physical interpretation of the delta-prime interaction is similar to the one of Seba ^H] 
but it does not require any renormalization of the coupling constant. 
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from the following momentum dependent interaction, 



V = c6{x) + 4XdJ{x)d^ + 2(7 + ir])dj{x) - 2(7 - ir])6{x)d, 



'X, 



(6) 



where c, A, 7 and rj are real coupling parameters. This interpretation is differ- 
ent from the standard one as a limit of conventional potentials which requires 
renormalizations of the coupling constants ^7j . However, it leads to a convenient 
parameterization (c. A, 7, 77) of the local interactions which has a natural physical 
interpretation and which is without constraints. Moreover, as shown in Section 
2, by simple formal computations the interaction in © is turned into standard 
boundary conditions which provide a mathematically rigorous formulation of the 
model. 

We can now give a more specific description of the many-body systems we 
consider. They are defined by the Hamiltonian 



with Xj G M the particle coordinates and Vjk local two-body interactions obtained 
from the momentum dependent potential in (jH)) by replacing x by the inter-particle 
distance Xj — Xk- In our discussion of special cases of this model we will use the 
notation (c, A,7,77) introduced above, e.g., {c, 0,0,7]) refers to the model defined 
in (P). It is important to note that the general model {c,\,'y,ri) describes iden- 
tical particles only for rj = j = (since the interaction terms dxS{x) ± 6{x)dx in 
(in)) change sign under particle exchang ). Thus, besides the 

delta-interaction model (c, 0, 0, 0), also the model (0, A, 0, 0) is of special interest. 
However, different from the former, the latter is not exactly solvable in the general 
case of distinguishable particles ^j. Still, the restriction of this model to fermions 
is interesting since it provides the only non-trivial exactly solvable fermion model 
with local interactions in ID (recall the the delta interaction is trivial for fermions; 
note that we only discuss models without internal degrees of freedom here). More- 
over, it has a natural physical interpretation as the non-relativistic limit of the 
massive Thirring model [12J. Our results in this paper imply that the model in 
is the only generalization of this fermion model to distinguishable particles which 
remains exactly solvable. 

As already mentioned, the very same questions studied in this paper was also 
studied by Albeverio et. al. [IHl,^ but their results differ from ours. The reason for 
this is that in ^Ul it is assumed that the model describes identical particles so that 
the approach in [5] applies, and this restricts their analysis (implicitly) to the two 
parameter family (c. A, 0, 0) of local interactions. Thus, while jTUI also finds that 

*We learnt about this work after finishing ours. 
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(7) 
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the model in Q is integrable, it concludes that the model is integrable only 
for the previously known case t] = 0.^ Indeed, one can check that the boundary 
conditions found to be integrable in (TU], Eqs. (15) and (16), correspond to the 
cases (c, 0,0,0) and (g, l/g,0,0) in our terminology.^ Thus, our main result in 
addition to ^U] is the extension of Yang's approach j3] to models of non-identical 
particles and thereby giving the first conclusive answer to the question posed in the 
second paragraph above. Moreover, we also give a novel physical interpretation 
of these models. 

The remainder of the paper is organized as follows. In Section 2 we discuss 
local interactions in ID, and in particular the relation of our interpretation to the 
standard one. In Section 3 we derive the consistency conditions on the coordinate 
Bethe Ansatz for the many-body systems with local interactions. Solving them we 
obtain the two families of systems defined in (P) and © (Section 3.2). Section 3.3 
contains the explicit solution of the latter model. In Section 3.4 we construct the 
recursion relations for the coefficients arising in the coordinate Bethe Ansatz for 
the model in which leads us to a somewhat unusual representation of the 
Yang-Baxter relations. We conclude with a few remarks on interesting related 
models not considered in this paper (Section 4). Some technical results related to 
our solution of the Yang-Baxter relation are deferred to two appendices. 

2 Local interactions in ID 

Interactions localized at points in ID have been studied extensively using the 
mathematical theory of defect indices; see [21 and references therein. From these 
studies it is well-known that the delta interaction is only one of many possible 
local interactions, and that a general such interaction can be characterized by 
four real coupling parameters. This can be formally understood as follows: for 
a ID Hamiltonian H = —d^ + V with an interaction V localized at x = 0, all 
eigenf unctions iIj{x) should be smooth everywhere except at x = 0, and {HiIj){x) = 
—ip"{x) for non-zero x. Requiring H to be self-adjoint leads to the following 
consistency condition. 



^In a system is referred to as integrable if it is exactly solvable by the coordinate Bethe 
Ansatz. 

^The third case obtained in ^H] corresponds to the limiting case lim^_»oo(c'7^, 0, 77, 0) in our 
terminology and thus violates the assumptions under which it was obtained. 





for arbitrary wave functions and if)^ or equivalently 




(9) 
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General boundary conditions are of the form 

^'{0+) = MnV^'(-0+)+Mi2Z^(-0+) 

^(0+) = M2iV''(-0+)+M22^(-0+) (10) 

(and similarly for 0, of course) and are thus parameterized by four complex param- 
eters Ujk which, upon imposing Q, are reduced to two complex, or equivalently, 
four real parameters. 

One prominent example of such consistent boundary conditions are 

^(0+) = ^(-0+) 
^^'(0+) -^'(-0+) = c^/'(0+), (11) 

which, as is well-known, corresponds to the delta interaction in Another 
example is what is usually referred to as delta-prime interaction, 

^'(0+) = ij'{-0+) 
V'(0+) -^^(-0+) = 4A^'(0+), (12) 

and which corresponds to the momentum dependent interaction in (0). 

As mentioned in the introduction, we use the following physical interpretation 
of the full four parameter family of local interactions, 

H = -dl + V, V = c6{x) +AXdJ{x)d^ + 2{-f + ir])dj{x) - 2{-f -ir])6{x)d^, (13) 

where c. A, 7 and 1] are real coupling constants. Once realized, this result is 
rather plausible: the operator V is obviously local and (formally) self-adjoint, 
and it contains four real parameters. Moreover, it is obviously the most general 
such interaction containing at most two derivatives, and it is plausible that higher 
derivatives than that cannot lead to consistent boundary conditions. It is also 
not so difficult to show that it indeed corresponds to the four parameter family of 
local interactions introduced in P^l^: formally, the eigenvalue equation of H in 
(HSl) is equal to 

- ij"{x) + c6{x)iIj{0) + AX6'{x)tlj'{0) + AX6{x)ij"{0) + 
+2(7 + ir/)5'(x)^(0) + 2(7 + iri)S{x)ip'{0) - 2(7 - iri)S{x)ip'{0) = Eip{x) (14) 

with ijj the eigenfunction, E the corresponding eigenvalue, and primes indicat- 
ing derivatives. The crucial point is that, due to the singular interaction, the 
eigenfunctions are in general not continuous at the singular point x = 0, and one 
therefore has to interpret the eigenfunction and its derivative at the singular point 
as the average of the left- and right limits, 

^(0) = l[^(o+) + ^(-0+)] 

^'(0) = ^[^'(0+)+^'(-0+)]. (15) 
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Using that ()14|) can be turned into boundary conditions of the form (fTUj) in the 
following way: integrating once from x = —0'^ to x = 0^ , and twice, first from 
X = —0+ to X > and then once more from x = —0~^ to x = 0"'' gives 

V''(0+) - ^/''(-0+) = ciIj{0) -2{-f -ir])^'{0) 
V^(0+) - V^(-0+) = 4AV^'(0) + 2(7 + ir7)^(0). (16) 

It is interesting to note that this parameterization of boundary conditions is iden- 
tical with the one proposed in jTB], Eq. (2.1). Using (jisp we can write these latter 
equations as in (fTUI). with 

u ^ = iu^ru^, u^ = ('^ ^i- ^ . ,] , (17) 

\u2i U22) ^ V T2A l^(7 + iry)y' ^ ^ 

where we introduce a convenient matrix notation.'' It is easily seen that the 
condition in ^ is equivalent to 

J=(° ~^^=U^JU, (18) 

which indeed is satisfied by the matrix U in ()17p . We note that this physical 
interpretation provides a convenient parameterization of the full four parameter 
family of local interactions. In particular, it is without constraints, in this respect 
distinguishing it from others used in the literature, e.g. jH HHI 

t/^e«(: ;■), ..[0..), M.«...Kwith3*-..^l. (19) 

Rephrasing our result in this latter parameterization is straightforward ^H], but 
its physical interpretation would be less clear. 



3 Coordinate Bethe Ansatz 

In this section we consider the model of particles interacting via the full four 
parameter family of local two-body interactions. It is defined by the Hamiltonian 
Hn in O with 

Vjk = 2cS{xj - Xk) + 2\{d^,^ - d^^) S{xj - Xk) {d^^ - d^,^) 

+ 2(7 + 17]) {d,^ - d,,) - 2(7 - i7])6{xj - Xk) {d,^ - , (20) 

^The special cases det([/+) — where U above is undefined are discussed in TB^. We only 
mention that the case = 7 = and A = 1/c corresponds to separated boundary conditions, 
V''(±0+) = c^(±0)/2. 
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which is the obvious particle generahzation of ©• Using 2dx -Xk = (^x ~dxf. it 
is straightforward to generahze the boundary conditions of the previous section to 
the A^-particle case. We thus deduce that the eigenfunctions ip of if at are defined 
by the free Schrodinger equation 

(^^d] + E^ij{xi,...,xr,) = for Xj ^ Xk, (21) 

and the following boundary conditions, 

{dx, - dx,) [ij\x^=x,+o+ - V'U,=Xfe-o+] = 2c^U^.=^, (22) 
-2(7 -i^) {dx, - dx,) 'ip\x,=x, 
i'\xj=x,+o+ - i'\xj=Xk-o+ = 2A {dx^ - dx^) ^\x,=x, + 2(7 + iri)^\xj=x^, (23) 

where j < k and where, similarly to the one particle case in the previous section, 
the eigenfunctions as well as their derivatives must be regularized at points of 
interaction follows 

i^\xj=Xk = 2 ['^\xj=Xk+0+ + i^\x,=Xk-0+] 

{dx, - dx,) Mx^^x, - \ {dx, - dx,) [^U,^..+o. + ^U,^..-o.] . (24) 

Equations ()21|) - (f^ provide a mathematically rigorous formulation of our general 
model. 



3.1 Two particle case 

Before continuing our discussion of the full particle model we consider in some 
detail the two particle case, N = 2, elucidating some of the properties that make 
this model rather special. 

The most important such property is that the Hamiltonian H2 is invariant 
under permutations of the coordinates xi and X2 only in the case 7 = r/ = 0, 
and that we therefore, in general, have a model of non-identical particles. To 
determine the implications this has on the eigenfunctions we start by considering 
the following Ansatz for the two particle scattering states, 

S^ih - A;2)e''=i^i+^^2X2^ X2<xi ^ ' 

which, upon substitution into boundary conditions (P^. result in the following 
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expressions for the two scattering amplitudes and S^, 

(7^ + _ 2iri + cA — l)u 



S+{u) 



iAn^ — (7^ + 77^ + cA + l)u — ic 



/ N iAw^ + 27M + ic 

'^^^''^ ~ iAn2 - (72 + ry2 + cA + l)^z - ic' ^ ^ 

where we have introduced u = ki — k2. For future convenience we introduce also 
the scattering amphtudes and obtained from S"^ and by reverting the 
sign of 7 and r], 

S*^ = S'_^|^^_^^^_>_^, Sj, = S^lj^^j^fj^^j^. (27) 

From this Ansatz we construct a second set of eigenfunctions of H2 by noting that 
the Hamiltonian H2, while not invariant under the exchange of the coordinates xi 
and X2, is invariant under the simultaneous exchange of the coordinates xi and 
X2 and replacement of 7 and rj by —7 and —t] respectively. This implies that also 
(I)\xi^x2,'y^-'y,v^-ri is an eigeufunctiou of H2, or equivalently, that (f)\x-i^^x2 is an 
eigenfunction of H2\^^-^,n^-r,- Note that (/> and (f)\xi^x2,'y^-'y,rj^-v linearly 
independent. Since the potential V12 = unless Xi = X2, every eigenfunction 
ip of H2 obeys the free Schrodinger equation {df + 9| + -E) ip{xi, X2) = in all 
regions without coinciding coordinates. Thus, in these regions they are linear 
combinations of the plane waves ^^^^^^+^'^2x2 ^j^^^ ^ik2Xi+ikix2 ^ This implies that 
and (f)\x-^^x2,'y^-'y,rj^-v provide, in the absence of bound states, a complete set of 
eigenfunctions of H2. 



3.2 Consistency conditions 

We now consider, for arbitrary A^, the model defined by in (jT)) and In 
particular, we determine for which values of the coupling constants (c. A, 7,77) its 
eigenfunctions iIj{x) can be obtained by the coordinate Bethe Ansatz 

iIj{x) = J2 ApiQW'"'""^ (28) 

in the wedge 

'■ ^Q(i) < ^Q{-2) < ■ ■ ■ < Xq(n) (29) 

with X = {xi, . . . ,xn) and kp ■ xq = Yl!j=i^P{j)^Q{j)-: Q E Sn The 

corresponding eigenvalue is obviously E = J2j'=i^'j- We recall that the validity 
of the Bethe Ansatz amounts to the model being quantum integrable in the sense 
that the most that can happen in any two-body scattering process is an exchange 
of particle momenta. 



9 



To take the boundary conditions into account we consider the boundary 
Xj = Xk for fixed j and k such that j < k. The last requirement is important since 
the particles in general are non-identical, as previously discussed. Furthermore, 
let Q be an element in Sn such that xg^i) = Xj and = Xk for some fixed 

i, implying that xj = — 0"*" is contained in the wedge Ag and xj = + 0"*" in 
AgTj, where Tj e Sn is the transposition interchanging i and i + 1. From this we 
deduce that the boundary conditions imply the following relations between the 
coefficients Ap{Q) of the coordinate Bethe Ansatz (|28|). 

i(fcp(i) - A;p(,+i)) [ApTAQT^) - Ap{QT,) - Ap{Q) + AptXQ)] = c[Ap{Q) 
+ AptM + ApiQT,) + ApT^QTi)] - (i7 + v) [ApTAQT^) - ApiQT,) 

+ Ap{Q)-AptM] 

Ap{QTi) + ApT^QTi) - Ap{Q) - AptM = '^K^Pi^ ' kp{^+l)) [ApTXQT^) 
- Ap{QT,) + Ap{Q) - AptXQ)] + (7 + i^) [Ap{Q) + AptXQ) 

+ Ap{QTi) + ApTXQT^)\. (30) 

These relations constitute a linear, homogeneous system of 2(A^— 1)N^? equations 
for the iV!^ unknowns Ap{Q), and the Bethe Ansatz is consistent if and only if 
this over-determined system of equations has A^! independent solutions where the 
Ap{I) can be chosen arbitrarily. In the following discussion we will show that this 
is the case in the two special cases corresponding to the models defined by (P) 
and (0). To do this we will not attempt to solve these rather complicated system 
of equations by brute-force, but rather use a somewhat indirect approach which 
essentially amounts to reducing the A^-particle case to a sequence of two particle 
problems. 

Lemma 3.1. Let Q e Sn and i e {1,...,A^}. If Q{i) < Q{i + 1) and the 

coordinate Bethe Ansatz is consistent, then 

AptAQ) = Sp {kp(i) - fcp{i+i)) Ap{Q) + 5*^ (A;p(i) - fcp(i+i)) Ap{QTi) 
AppXQTi) = Sp {kp{i) - kp(i+i))Ap{QTi) + S^{kp(i) - fcp{i+i)) Ap{Q), (31) 

where Ti is the transposition interchanging i and i + 1. 

Proof. We first consider the case N = 2. Let ip be an arbitrary eigenfunction of 
H2. Recall that the eigenfunctions and 01x1^x2,7^-7.'?^-'? -^2 constitute, in 
the absence of bound states, a complete basis for the eigenspace of H2. It follows 
that is a. linear combination, 

i> = ai(f) + a20|xi^X2,7,r,^-7 -r, 

r aie'^^i^i+^^^x^^ j^(^^_^2)e''=2xi+ifcix2^ xi < X2 . . 

\ /2('tl - 't2)e'''2^2+ifcixi ^^^gifciX2+ifc2Xi^ x^<xi ' ^ ' 
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for some complex constants ai and 02; we have introduced the functions 

fiih - k2) = S;^{ki - k2)ai + S:^{ki - /c2)a2 
f2{ki - k2) = S]^{ki - ^2)02 + S:^{ki - k2)ai. (33) 

Relabehng the constants ai and 02 as well as the functions /i and /2 as follows, 
ai = a2 = Aj(Ti), At^{I) = f\ and At^{Ti) = f2, we arrive at the two- 

particle Bethe Ansatz with coefficients Ap{Q) given by 

^T,(/) = S^{k, - k2)Aj{I) + ST{k^ - k2)Aj{T,) 
AtATi) = S^ik, - k2)AjiT,) + S+ik, - k2)AiiI). (34) 

To extend this result to arbitrary values of we observe that the boundary 
conditions (|22j) become identical with the boundary conditions for the case N = 2 
if we substitute Xj — > Xi and Xk — > X2, for all j < k. For each P G S^, let 
Aq be a wedge such that xqi = Xj and a;Q(i+i) = Xk- From the equivalence of 
boundary conditions, in the sense stated above, follows that the relations between 
the coefficients Ap{Q), Ap{QTi), AptXQ) and ApT^QTi) are obtained from (jH^ 
by the substitutions Xi ^Qii) and X2 — >■ XQ^i+i^ as well as ki kpi^i) and 
^2 kp(i+iy This yields the relations in (jHH). □ 

It is important to note that there is a possible inconsistency in the coordinate 
Bethe Ansatz arising from the fact that the representation of an element in Sn is 
not unique. However, any two representations can be converted into each other 
by using the defining relations of Sn, 

Ti rri -| rri rri rri rri rri rri 

i-Li = J-, J-i-Li+l-Li = -Li+lJ-i-Li+l 

T,Tj = TjT, for|z-j|>l. (35) 
Thus no inconsistency can arise provided that 

APT^T, (Q) = Ap{Q), ApT^T^^^T. (Q) = ApT^^^T.T,^^ (Q) 

ApT^T, (Q) = Apt^tM for - jI > 1 (36) 

for all P,QeSN. 

To determine when these conditions are valid we follow the approach of A. B. 
and AI. B. Zamolodchikov [201 and make use of the algebraic structure which has 
come to be known as the Zamolodchikov algebra. We start by briefly recalling the 
construction of the Zamolodchikov algebra: to each particle of type A and with 
momenta k the symbol A{k) is associated, and the two particle scattering theory 
is encoded in the commutation relation 

A{k,)B{k2) = Si''{ki2)A{k2)B{ki) + S^''{ku)B{k2)A{k,), (37) 



11 



where ku = ki — k2 and and S*^ are the two particle scattering amph- 
tudes. The full N particle scattering theory is then obtained by factorizing each 
scattering event into a product of two particle events. As observed in iden- 
tifying each product A{ki)B{k2)C{k3) . . . with a particular coefficient Ap{Q) of 
the coordinate Bethe Ansatz the consistency conditions ()36|) are equivalent to 
requiring the Zamolodchikov algebra to be consistent as well as associative. In 
our particular case, since the particles are distinguished only by their relative or- 
dering, it is sufficient to consider the case of three particles. Therefore, let the 
product A{ki)B{k2)C{k3) correspond to a particular coefficient Ap{Q), for fixed 
P,Q & S3. A straightforward but somewhat tedious computation then shows that 
the Zamolodchikov algebra, in our case, is consistent as well as associative if and 
only if the two-particle scattering amplitudes and obey the following so 
called Factorization equations, 
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(47) 


si 


[v)Sj^{u + v)Sj^{u) 


+ 5^ 


iv)St 


{u + v)S^{u) 


= Sr 


u)Sq^{u + v)S^{v) 


(48) 


Sr{'")Sr{u + v)S^{u) 




iv)SR 


{u + v)S-^{u) 


= Sr 


u)Sj^{u + v)Sj^{v) 


(49) 


Sr 


['^)Sr{u + 'v)S^{u] 






{u + v)S'^{u) 


= s^ 


u)S^{u + v)S]^{y) 


(50) 



for all real u and v. 

Upon substituting the two-particle scattering amplitudes and into the 
Factorization equations above a straightforward but somewhat tedious calculation 
shows that (jSHl)-(IlH) as well as ^^-^^ are fulfilled for all values of (c, A,7,r7), 
while (jl2I) and (Iini)-(ISni) holds true if and only if 

7[A(3c + \{u^ + UV + v"^)) + 47^] = 
[-1 + cA + 7^ + ri{ri - 2i)] [A(3c + \{u^ + uv + v"^)) + 47^] = 
[-1 + cA + 7^ + r]{r] - 2i)] [A(3c + X{u^ + uv + v'^)) + 47^] = (51) 

for all real u and v. Using the fact that the first condition requires 7 = and 
that the last two are related by complex conjugation, we can reduce them to the 
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rather simple form 



7 = 

A[-l + cA + r/(r7 - 2i)] = 0, (52) 

which obviously has two families of solutions, each of which defines systems exactly 
solvable by the coordinate Bethe Ansatz. We thus arrive at the main result of 
this paper. 

Theorem 3.2. The coordinate Bethe Ansatz for the eigenfunctions of the 
Hamiltonian is consistent if and only if 

A = 7 = (53) 

or 

A = 1/c and 7 = = 0. (54) 

Note that these two sets of conditions on the coupling constants (c, A, 7, rj) 
correspond to the two Hamiltonians defined in and It is interesting to 
note that, in both cases, 5*^ and 5*^ only have a single pole, and 5^ = 5*^ even 
in the first case. 



3.3 Explicit results for the case (c, 1/c, 0,0) 

This particular case correspond to the Hamiltonian H^'^^ in Eq. (j2)). We note 
that this Hamiltonian is invariant under permutations of the coordinates Xj, and 
one can therefore assume a particular exchange statistics, which determines the 
eigenfunction in all wedges once it is known in one of them. Furthermore, using the 

fact that 5'^|a=i/c,7=j7=o = "S'j. |a=i/c,7=?7=o = and S^\x=i/c^^=ri=o = 5'^|a=i/c,7=77=0) 
the relations in (j31|) reduce to the rather simple form 



ApTiiQ) = 5'^|A=i/c,7=r)=o {kp{i) — A;p(j+i)) Ap{Q) 
_ i (fcp(i) - A;p(i+i))^ + ic^ 



-Ap{Q) 



i {kp(^ij - fcp(i+i)) + 2 {kp(i) - fcp(i+i)) - ic' 

= 't"~'T'*i:' MQ) (55) 

1 (fcp(i) - fcp(i+i) ) + c 

for all Q G 5'iv- This implies that 

ApTi = Yi [kpii) - kp{i+i)) Ap, (56) 
where we have introduced the function 

y.(u) ^ |±£. (57) 
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and interpret Ap as a vector with A^! elements Ap{Q). It is important to note 
that the functions Yi, for arbitrary exchange statistics, are identical to the cor- 
responding functions appearing when restricting the delta-interaction model to 
bosons (see e.g. The eigenfunctions of the latter model are well-known (see 

e.g. Section I.l in 



ij = C 



.N>j>k>l 



det [exp(iA;„a;„)], 

l<m,n<N 



Xi < X2 < 



(5^ 



where C is a normalization constant. We conclude that this formula also gives 
the eigenfunctions of the Hamiltonian in (j2)) in the identity wedge A/, for arbi- 
trary exchange statistics (characterized by some Young tableaux). Using standard 
arguments from group theory (see e.g. |M|); these eigenfunctions can be straight- 
forwardly extended to all other wedges. 



3.4 Recursion relations for the coefficients A 



p 



In order to provide a machinery for computing the eigenfunctions explicitly we 
now proceed to derive a set of recursion relations for the coefficients of the 
coordinate Bethe Ansatz. Our starting point is the following fact: defining 



(59) 



one can write 



Ap{QR) = J2 iR)Q,Q'MQ') = mpKQ), 



(60) 



where the first equality is a trivial consequence of the definition, and in the second 
we interpret {R)q^qi as elements of an n x n matrix R with n = N\ the rank |S'7v| 
of Sn- These matrices obviously define a representation R ^ R of Sn acting 
on the coefficients Ap{Q). It is worth noting that this is identical with the so 
called (right) regular representation of S^. Using this fact we can rewrite the two 
relations in (j31|) as follows, 

AptAQ) = S;^ (kpii) - kp(^i+i{) + {kp(^i) - kp^i+i))fi Ap{Q) 
AprAQTi) = S]^ {kp(i) - kp(i+i)) + {kp(i) - kp(i+i)) fi Ap{QTi), (61) 

where Q is required to fulfill the condition Q{i) < Q{i + 1). To organize the 
elements of the vector we have to introduce an ordering on the set of permu- 
tations. We find it convenient to use the following 
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Definition 3.3. Associate to any two permutations Q, Q' G Sn the sequence 
Qi := Q{N — i + 1) — Q'{N — i + 1), where i = 1,2, . . . , N . If the first non-zero 
number in the sequence {aj} is positive, Q is said to be larger than Q' , denoted as 
Q > Q'. 

The idea for the remainder is the following: order the coefficients Ap{Q) ac- 
cording to the ordering just defined (largest permutation first) into the vector Ap, 
and use the recursion relations in for the individual coefficients Ap{Q) to 
write 

Apt, = Yi {kp(i) - A;p(i+i)) Ap, (62) 
where we have introduced the matrix 

Yi{u) = S^m) + S^(M)fi (63) 

in which S^j and are diagonal A^! x A^! matrices with entries and 
respectively. To determine the distribution of the scattering amplitudes 5*^ and 

among the diagonal elements of the matrices S\ and we start by deducing 
a natural decomposition of an arbitrary permutation into a product of elementary 
transpositions Tj. This decomposition will provide us with enough information to 
determine the structure of and explicitly. 

We pause to introduce some notation and conventions to be used in the remain- 
der of the discussion. By abuse of notation, we identify a permutation Q E 
with Q' G S]\f+n, n > 0, a 

Q{i)=Q'{{} for alH = l,2,...,Ar 
Q'{i)=i ioT alii = N + 1,N + 2,. . . ,n. (64) 

E.g., (231) G 5*3 and (23145) G will be identified.® Furthermore, we denote 
as the kth permutation of Sm, ordered according to the ordering defined in 
Definition HO E.g. for ^3 this implies that Qf = (123), Ql = (213), Q| = (132), 
Ql = (312), Ql = (231) and Ql = (321). 

Using the fact that every integer j such that 1 < j < m\ can be uniquely 
written in the form j = n{m — 1)1 + k, where < n < m — 1 and 1 < A; < (m — 1)!, 
we prove 

Lemma 3.4. Let k, m and n be integers such that m>l,0<n<m — 1 and 
1 < ^ < (m — 1)!. Then 

m— 1 

Qn(m-l)!+fc ^ n ^' ^^^^ 

i=m—n 

^We recall that a permutation Q = (ijk) is defined such that Q{1) — i, Q{2) = j and 
Q(3) = k, etc. 
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where 



Wt, = \ Y'^'---^"' . . (66) 

-LA ntheriinse ^ ' 



otherwise 



Proof. It clearly follows from Definition 13 . 31 and the fact 15*^1 = rn\ that Q'^(^^_iy_^_i 
is obtained by cyclically permuting the last m—n elements of the identity permuta- 
tion in Sm such that Q^(^^_iy_^_i{m) = m—n. Decomposing this cychc permutation 
into a product of elementary transpositions Tj, we obtain 

Qn{m-iy.+l ~ Tm-n ■ ■ ■ (67) 

We conclude the proof of the Lemma by observing that 

Q^{m-iy.+k = Q^im-iy.+iQT ^- (68) 

□ 

Corollary 3.5. Let k, m and n be integers such that m>l,0<n<m and 
1 < k <m\. Then 

Qn:^+M) > Qn:^+M + 1), l<k<n{m- 1)! 
e+fc(^)<^?™!+fc(^ + l). n{m-l)\<k<m\. (69) 

Proof. From Lemma 13.41 clearly follows that Q™^|^^(m + 1) = rn + 1 — and 
consequently that 



e+i(j)<Q™^!^+i(^ + l), l<J<m 



' nm!+l 

Qn^w(j)>QI^w(^ + l), m-n<j<m. (70) 

In other words, for each fixed n there exists n distinct integers I < m + 1 such 
that / > Q'^^i^^^^m + 1). This, together with the facts that Q™+/+;, = Q'^^^+^Q'^ 
and that Q^{m) = m, for all k < {m — 1)!, imply the statement. □ 

Note that Lemma 13.41 provides a recursive procedure for decomposing an ar- 
bitrary permutation into a product of elementary transpositions Tj. 

Using Lemma 13.41 and Corollary 13.51 we now determine the structure of the 
matrices S^j and S^. Corollarv 13 . 51 and the recursion relations ()61|) together imply 
that 

_ j I <j -nil <n{i~l)\ 
^^^j \ S+ n{i-l)\ < J -nil <i\ ^ ' 



n 



and 



_ } 1 < J - < ^(^ - 1)! /79N 

1 St n{t-iy.<j-m\<t\ ' ^ ' 
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where n is required to fulfill the relation j = nil + k, for some 1 < k < i\. This 
determines only the first {i + 1)! diagonal elements of the A^! x A^! matrices 
and S^. However, it immediately follows from Lemma 13.41 that this structure is 
periodic with period {i + 1)!, i.e., 

(ST)j+(i+i)!j+(i+l)! = (St)jj5 (73) 

thus determining all A^!^ elements of and S^. We recall that for our case, 

Q±( \ (-^^ ± 2ir] + l)u ± -ic 

[r]^ + l)u + ic [1]^ + l)u + ic 

but we give the construction for the more general case where 5*^ are different since 
we hope that this will be useful for other models. 

In deriving the recursion relations ()62|) it is important to note that there is a 
possible inconsistency in the coordinate Bethe Ansatz, arising from the fact that 
the elementary transpositions Tj obey the defining relations ()35|) of the permuta- 
tion group Sn, as discussed in Section 3.2. Thus the recursion relations (jH^ are 
consistent if and only if 

ApTiT, = Ap, ApT^T,+iT, = ^pr,+ ir,T,+i, 

ApT^Tj = ApT^T, for |z - j| > 1 (75) 

for all P G Sn- Using the recursion relations (j62j) one find that that these latter 
conditions holds true if and only if the following so called Yang-Baxter relations 
are fulfilled, 

Y,{-u)Y,{u) = /, Y,(t;)Y,+i(n + v)Y,{u) = Y,+i(n)Y,(n + v)Y.,^^{v), 

Y,{u)Y,{v) = Y^{v)Y,{u) for \i-]\>l (76) 

for all real u and v. However, since we already in Section 3.2 determined for which 
values of the coupling parameters (c, A,7,r7) the relations (f7H|) holds true we can, 
in view of Theorem 13.21 immediately conclude that we have proven 

Proposition 3.6. // and only if the coupling parameters {c,X,'y,T]) satisfy the 
conditions in l\5!^} or |^]), then the Yang-Baxter relations (|y6| ) are fulfilled, and 

Ap = yp{k)Ai, (77) 

where yp{k) is a product of the matrices Yj(A;p(j) — A;p(j+i)) obtained by repeatedly 
using ^U^. 
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Remark 3.7. In the special case 7 = r/ = the Hamiltonian is invariant 
under permutations of the coordinates, and in consequence the scattering ampli- 
tude 5*^1^=^=0 = >S'^l7=r;=o =• Sji and similarly 5*^1 ^=^=0 = S^\.y=jj=o =: St- A 
particularly interesting special case is the delta interaction model A = 7 = = 0, 
in which case the matrices Yj take the well known form 

V f w iufi + ci 

Y,{u)\x=^=^=o = —. , (78) 

m — c 

originally introduced by Yang f^. 

To illustrate the general discussion of this section we provide in Appendix A 
a more explicit account of the three particle case N = 3. In Appendix B we 
outline a simple, direct proof of the Yang-Baxter equations for the cases stated in 
Proposition 3.6. 



4 Concluding remarks 

In the introduction we argued that the Hamiltonian defined in (jZj) and ()2()j] is 
the N particle generalization of the most general Hamiltonian H = —d^ + V with 
local interaction V. However, not all A^-body Hamiltonians with local two-body 
interactions can be obtained in this way. A particular example of such a model is 
the quantum version of the derivative nonlinear Schrodinger equation [21] defined 
by the Hamiltonian 

N 

j=l i<k 

where is a real coupling parameter. This model does not fall in the class studied 
in this paper since the interaction depends also on the sum of the momenta of 
the particles involved in the interaction, and it is therefore not Galilean invariant. 
However, this Hamiltonian is nevertheless interesting. It describes identical parti- 
cles, and it was shown by Gutkin [7j that it is exactly solvable by the coordinate 
Bethe Ansatz only in the case of bosons or fermions. 

All possible local two-body interactions are formally defined by the 2-body 
Hamiltonian H2 = ~ + ^12, where iWi2ip){xi,X2) = for two-particle 
wave functions ip except in regions of coinciding coordinates, Xi = X2- We believe 
that there is a nine parameter family of such two-body interactions which one 
should be able to find using general methods discussed in [21 CH] • It would be 
interesting to know if there are additional such distinguishable particle models 
exactly solvable by the coordinate Bethe Ansatz. 
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As discussed in jTH] Remark 2.4, there is a one parameter family of local inter- 
actions which are unitarily equivalent to the non-interacting case. This suggests 
that there should be also a one parameter extension of local interactions which 
are unitarily equivalent to the delta-interaction, and this might provide a sim- 
ple explanation of our results. At closer inspection we found that this is indeed 
the case: if and only if iIj{x) obeys the boundary conditions in ()16p defining the 
interaction {c,0,ri,0), then 

^(x) = e-^"®("V(a;), e'" = (80) 

1 — IT] 

with the Heaviside function B obeys the boundary conditions defining the interac- 
tion (c/(l -|- 7]"^), 0, 0, 0). In a similar manner, the model formally defined in is 
unitarily equivalent to the delta-gas with coupling c = c/{l + rf), and the unitary 
operator U intertwining the two models is given by 

(?7V')(xi, ...,xn)= e-^"^^<^ ...,xn) (81) 

with a as above. ^ Moreover, it is interesting to note that ()74|1 can be written as 

S^{u) = e^'%{u), S^{u) = a{u) (82) 

where Yi{u) = a{u) + h{u)Ti gives the rational solution of the Yang-Baxter rela- 
tions, suggesting that our solution Yi{u) of the Yang-Baxter relations is unitary 
equivalent to the latter. We should mention that this very argument was used 
already in ^1] to deduce the integrability of the model in However, while 

this is a simple alternative argument showing integrability of the model in (^, 
only our argument allows to conclusively determine all integrable cases. 

To conclude, in this paper we found an example of a ID quantum many-body 
system of non-identical particles which is integrable, and it led us to an interesting 
class of solutions of the Yang-Baxter relations. While the solution which we found 
is closely related to the well-known rational one, there might be others which are 
truly different from known ones, and our results provide a convenient starting 
point for finding them. It thus is tempting to speculate that our model is only a 
first example in a novel interesting class of integrable systems. 
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Appendix A. Recursion relations for the three 
particle case 

To illustrate the abstract discussion of Section 3.4 we here provide explicit formulas 
for the three particle case. We start by noting that the vector Ap, in the ordering 
defined in Section 3.3, has the following structure, 

A^p= (Ap(123), Ap(213), Ap(132), Ap(312), Ap(231), Ap(321)), (Al) 

where A^p denotes the transpose of Ap. Straightforward computations (or alterna- 
tively Lemma show that we can write this using the elementary transpositions 
Ti and T2 as follows, 

A'p={Ap{I), Ap{T,), Ap{T,), Ap{T,T,), Ap{T,T,), ApiT.T^T,)). (A2) 

Directly applying the recursion relations ()6H) we find that the matrices Sp and 
S^, in the particular ordering we have chosen, take the following form. 



(S^p 
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(A3) 



Using T1T2T1 = T2T1T2 and (jUTj) we deduce that 
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^0 














S}) 



(A4) 

It is straightforward to verify that this structure is indeed reproduced also by 
the general discussion presented in Section 3.3, in particular Eqs. (f7T |) - (f7H|) . To 
explicitly construct the matrices Yj, for i = 1,2, and thereby also the recursion 
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relations (jU^ . there remains only to determine the explicit form of Ti and T2. 
Using (j^UI) we find that 



/o 


1 











o\ 


1 


























1 














1 


























1 













1 


V 



To 



/O 1 0\ 

1 

1 
1 
1 

\0 1 0/ 



(A5) 



Inserting the above matrices into the consistency conditions ()76|) a straightforward 
but somewhat tedious calculation shows that they are equivalent to Factorization 
equations (jHS|) -(f ^ . This reflects the fact that, for this model, it is sufficient to 
consider the three-particle case in order to find those cases exactly solvable by the 
coordinate Bethe Ansatz. 



Appendix B. Direct proof of Yang-Baxter equa- 
tions 

In this Appendix we outline an alternative, direct proof of the validity of the 
Yang-Baxter relations in (f7n|) for the operators Yi{u) defined in Eqs. (fHHj) and 
(p3|) -P7 j) and the cases stated in Theorem 3.2. 

We first note that the formulas given in Appendix A allow a simple verification 
of these relations in the simplest non-trivial case = 3.^^ The general result, for 
arbitrary N, then follows from the following 

Lemma B.l. // the operators Yi{u) satisfy the Yang- Baxter relations for N = 3 
they satisfy them for any N > 2. 

This is true since all identities in fl76|) can be brought to a block diagonal form, 
where the relations in each block are unitarily equivalent to the corresponding ones 
for the case A^ = 3. 

To be more specific, we consider first the second identity in (j76|) . involving 
three matrices Y on each side, for some fixed N > 2 and i < N. It follows 
from (jU^ that Yi{u) acting on Ap mixes only elements Ap{Q) and Ap{QTi) 
(since S^j^^ are diagonal matrices). Thus, for any Q G Sn, both sides of the sec- 
ond set of identities in (fTBj) mix only the elements Ap{Q), Ap{QTi), Ap{QTi+i), 
Ap{QTi+iTi) and Ap{QTiTi^i), i.e., it is possible to reorder the elements of the 
vector so that one can group them into blocks of six elements such that the 

^^The verification can easily be performed by using a symbolic programming language like 
MAPLE or MATHEMATICA. 
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identities decompose into blocks of 6 x 6 matrices. Moreover, let Q' be the 
largest element, when using the ordering defined in Definition 13.31 in the set 
{Qi QTi, QTj+i, QTi^iTi, QTiTi^i, QjTj+iTj}, then 

Q' > Q'T, > Q'Ti+i > Q'T,+iT, > Q%T,+^ > Q%Ti+,T,. (Bl) 

This implies that the matrices Yj(u) and Yj_|_i(u) restricted to the vector 

{Ap{Q'), Ap{Q%), Ap{Q%+,), Ap{Q%,+,T,), Ap{Q%T,+{), Ap{Q%T,+{r,))' 

(B2) 

are identical with the 6x6 matrices Yi(m) and Y2(m) for = 3 respectively 
(these latter matrices are explicitly given in Appendix A). This proves the lemma 
for the second set of identities in (|7BJ. The verification of the first and last set of 
identities follows from similar arguments (they reduce to blocks of 2 x 2 and 4x4 
matrix identities which are easy to verify). 
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